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Abstract 

Domain walls of a discrete model of an anisotropic ferromagnet are studied. They 
can be described by sequences of two reversible mappings. Competition between the 
length scale of spatial structures and the lattice constant leads to a rich diversity of 
domain wall solutions related in a bifurcation scenario. 
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Introduction 



The complexity of spatial structures of discrete oscillator chains raises questions of consid- 
erable interest [l]-[8]. In this paper domain wall solutions of an one-dimensional discrete 
model of a ferromagnet are contrasted to the well-known Bloch wall solution of the corre- 
sponding continuous system. As is known 0], stationary planar solutions of 



Sx{J^S + aS,e,) (1) 



are represented by the equation (p" = —j cos(/>sin(^ of a point mass in a periodic potential; 
beside spatially periodic solutions, the kink solution 

sin (j) = ibtghw (2) 

V " 



is obtained. However, if the length y ^ of this kink is not large enough, the continuum 
approximation [|10| fails and one expects competition between the length scale of a spatial 
structure and the lattice constant. Several papers have studied the corresponding classical 
discrete system Q, Q. 

This paper is devoted to study of stationary localized structures connecting two homo- 
geneous domains in a model of a ferromagnetic chain with one kind of atom; quantum 
corrections @ are neglected. 

Stationary solutions of oscillator chains with appropriate next neighbour coupling can be 
described by mappings; e.g. the Prenkel-Kontorova system is associated to the standard 
map 0,||5|. In contrast, complexity of spin chains arises from the fact that the coordinates 
of one spin oscillator are not given unambiguously by the position of its two predecessors. 
For planar solutions two positions of the following spin are possible. Consequently only a 
small part of the stationary states can be reached by a single map. At every single step 
one can decide which one of two possible maps is to be applied, so that the system may be 
considered an explicitly time-dependent dynamical system. Thus, according to symbolic 
dynamics, a subdivision of the solutions with their defining symbol sequence suggests itself. 
Solutions which are characterized by the same inversion symmetry as the Bloch wall are 
related to symmetry lines of the mappings. 



The anisotropic spin chain 

We discuss stationary planar solutions of the model of a spin chain Q , , Q with isotropic 
ferromagnetic coupling and an uniaxial anisotropic term 

Sn = SnX (J(5„+i + Sn-l) + aSnzSz) (3) 

J is the isotropic Heisenberg coupling constant, a an anisotropy constant; these constants 
are real. The modulus of the spin is conserved, we set |S'„| = 1. S„ is fixed if the local 
effective field /le//™ = J{Sn-i + Sn+i) + aSnz^z is parallel to Sn- The position of Sn+i 
can be computed from the positions of Sn and Sn^i- This determination is not unique; 
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in general several directions Sn+i yield parallel effective fields heffn- In the same way, 
possible directions of Sn-2 can be computed. 

The system is rotational symmetric with regard to the z-axis. Introducing spherical coor- 
dinates 

Sn = (cos 9n, sin 9n sin Cpn, sin 9n cos (j)n) (4) 

we investigate solutions with spins perpendicular to the x-axis (9 = In this case, fixed 
point conditions are given by the system of equations 

sin(</)„+i - (j)n) = sm{(j)n - (pn-i) + jj sin(20„) (5) 

Introducing the second coordinate In = sm{(j)n — (pn-i) the fixed point condition can be 
satisfied equally by 

0n+i = 4'n + arcsin In+l 
Li : In+i = In + k sin 20„ (6) 
k 



and by 



2J 

(j)n+i = (/)n + TT - arcsin 
L2 : I„+i = In + k sin 2(?!)„ (7) 

U — JL. 
— 2J 

Stationary solutions are defined by their initial conditions and a sequence of Li's and L2S, a 
property reminiscent of symbolic dynamics. In general, there are several solutions with the 
same sequence but with different initial conditions. In this sense, our system is equivalent 
to a lowdimensional dynamical system with an explicit time dependence which causes the 
sequence of mappings Li,L2. 

In must not exceed 1. For this reason, not all initial conditions lead to a stationary solution. 
The domain of definition is in general a complicated subset of the cylinder [0, 1] x [0, 2it]. 
Each mapping commutes with 

and with 

M- (9) 

Li and L2 are reversible mappings i.e. they can be decomposed into involutions 
L = HoG with = = Id (identity mapping). It follows that = GoH. Employing 
M~ and they can each be decomposed into involutions in four different ways: 

Li = o G°+ = Hl"^ o = H^- o G^- = H^- o G^- (10) 

L2 = Hl+ o G°+ = o G^+ = H^- o G°- = Hl^- o G^- (11) 
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We draw up the involutions and symmetry lines Fix{G),Fix{H): 



in vnln tinn 


SVlTlTTlPtrV llTlPS 


„o+ (f)' = —4> + arcsin J 
^1 '■ p = j 




^0+ . 4>' = -<P 

' 7' — / + A- sin 2(6 


4) = Q 


(p' = TT — (f) + arcsin 7 
^1 ■ p = i 


^ TT+arcsin/ 


. 0' = TT - (/> 

■ /' — 7 + A- sin 2fl!) 




rjQ- 0' = ^ — arcsin 7 
1 • 7' = -7 


7 = 


■ 7^' — — 7^ — A- sin 2r6 


T fe sin 2</) 
~ 2 


„7r- 4>' = 4> ~ arcsin 7 + tt 
1 ■ 7' = -7 





' I' = -I - fcsin2(i 


— 


0+ = — o!) + vr — arcsin I 


,i TT — arcsin/ 
2 


0' = — TT + arcsin 7 
^2 : /' = _/ 




„7r+ 1^' = — ^ + 27r — arcsin 7 

^2 • = / 




= + arcsin 7 
^2 : /' = _/ 


7 = 



The fact that Li and L2 have the same involutions G^/'^^ has a consequence of importance 
in the further argumentation. We assume U to be any finite sequence of Li's and L2's and V 
the corresponding sequence with reversed order of indices, for instance U = L2 o L2 o Li o L2 
and V = L2 o Li o L2 o L2. It follows that 

U-^ = GoV oG (13) 

(we skip the indices of G and 77 if they are unimportant). 



Continuous families of domain walls 

In this chapter symmetry properties of domain wall structures are studied. In our context 
we understand a domain as a homogeneous or spatially periodic structure. A domain can 
be defined as a fixed point of a finite sequence D = Li" o o Li^ ... . For example, 
(7 = 0, = 0) is a fixed point oi D = Li. This corresponds to the ferromagnetic state with 
parallel spins. (7 = 0, = 0) is also a fixed point of 7) = L2^. In this case it corresponds 
to the antiferromagnetic state of antiparallel spins. 
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Domain walls are stationary structures connecting two homogeneous domains. Let Di, 
D2 denote two finite sequences of Li and L2 which characterize the domains. [Idi,4'Di)-, 
{Id2,4'D2) are fixed points of Di and D2 respectively. A domain boundary solution con- 
necting the two saddle points {Idi,4'Di), {Id2,4'D2) may contain a sequence of mappings 
different from Di and D2, so a domain boundary solution is defined as a point set deter- 
mined by the sequence Di o B o D2 ■ Here Di and D2 denote infinite sequences of Di and 
D2 while B is an arbitrary finite sequence of Li and -L2. Such a domain boundary solution 
is obtained if one considers a point sequence emanating from a fixed point of Di on its 
unstable manifold W}^_^ . After the application of a certain sequence B of the mappings Li 
and L2 the point sequence must follow again the stable manifold Wl)^ towards a fixed point 
of D2; consequently domain walls are described as heteroclinic orbits of the saddle points. 
In correspondence with the solution of the continuous system, we are particularly interested 
in inversion symmetrical solutions. Such solutions can only be obtained if Di and D2 are 
defined as sequences with reversed order of indices, for example Di = Li o {L2)^,D2 = 
(L2)^ o Li or Di = D2 = L2. B can be an arbitrary finite symmetric sequence, for instance 
(1/2)^ or L2^ o Li o o Li o ^2^. (Iq, 4>o) is a point on the unstable manifold Wjj^ of the 
saddle point (IdiAdi) = L{°^{Io,(j)o). 

SJ,^ (Jo, 0o) = {...(/o, 0o), ihAi), •••} C W}^^ (14) 

is the set of all images Di{Io, (pQ),n G Z of {Io,(po). Applying an involution G to this set 
and using D^^ = G o Di o G (jl|), G^ = Id 

% = G{Sl^ (/o, M) = {...G{IoAo), G(/i, 0i), ...} (15) 

turns out to be the set of iterations {D2^^)"'{G{Io,(f>o)), which is an orbit on the stable 
manifold W^,^ of the saddle point G{Idi 1 0Di )• So the involution G transforms the unstable 
manifold of a fixed point into the stable manifold of another fixed point. 
The sets Sj^_^ and Sf^^ represent a heteroclinic orbit if they are identical. This is true if 
a common point of Sjj_^ and Sjy^ exists, e.g. a fixed point of an involution. S"]^^ or ^l^^ 
represent an intersection point set of Wj^_^ and W^,^ . This is the most simple case without 
a central sequence B different from Di or D2- 

For a domain wall associated to a sequence B, a point {Ib^^I^b) G must exist which 
is mapped by B onto W^)^. So we must find an intersection point of BiW"^^) and W^)^ in 
order to obtain a domain wall solution. B must cause a jump from a point of to the 
corresponding point of S'l)^; i.e. for one point (Is, (j)B) oi the conditions 

B{lB,(t^B) = G{lBAB) (16) 

must be fulfilled. We will show that sequences Di o B o D2 lead to structurally stable 
domain wall solutions. There are three different forms of the central symmetric sequence 
B: 

B = UoV (17) 

B = UoLioV (18) 

B = UoL2oV (19) 

where U and V have inverse order of indices (|T3|). We will show that any intersection point 
of ViW) or (|1|) Li o ViW) or (|l|) L2 o y(VF") and a symmetry line of (0) G, 



5 



Hi or (^) H2 yields an inversion symmetric solution. So to speak, we just have to find an 
intersection point of two lines on a plane; no additional requirement must be satisfied for 
reversibility. 

For a symmetric solution we have to find a point (/_b,0b) of VF" with the property (16) 
We focus the case ( [l7|) B = U oV . It will be shown that condition ( |l6|) follows from the 
existence of an intersection point (Ic", (pc) = ^(-^_B> <Pb) of V{W']j^) and a symmetry line of 
G^/'^^. Clearly such a point is invariant under the G concerned. With the property (p^ 
= G o U o G we have 

(Ib, (/>b) = V-\lc, (kc) = GoUo G{Ic, (t>c) (20) 

Applying G and using G^ = Id, G{Ic, (pc) = {Ic, 4>c) we obtain 

G(/b, Pb) = U{Ic, Pc) = Uo V{Ib, Pb) = B{Ib, 0b) (21) 

which had to be proved. 

For the case (^) B = U o LioV we have to prove 

U oLioV{lB,(t>B) = G'^'^^{lBAB) (22) 

where {IbiPb) is a point on W^_^. Presupposing an intersection point (IcPc) = Li o 
V{Ib, 4>b) of a symmetry line of H^^^^ and the line Li o V{W'^), we see that 

{Ic, 0c) = ^i(/c, 0c) = ^1 o ^i(/c-i, 0c~i) = G(/c-i, 0c-i) (23) 
and consequently 

{Ib, Pb) = V-\lc-i, 0c-i) = GoUo G(/c-i, 0c-i) (24) 
Applying G we obtain 

G{Ib, (Pb) = U{Ic, Pc) = UoLio V{Ib, Pb) (25) 



This can be carried out in the same way for (|T9|) B = U 0L20V, in which case we presuppose 
an intersection point of a symmetry line of H2^'^^ and the line L2 o V{Wj^_^). We conclude 
that symmetric solutions constitute a continuous family with respect to the parameter k. 
This proof does not imply that all solutions are symmetric. Intersection points of B{Wj^^) 
and Wl)^ may exist which do not correspond to any symmetry line. 



Bifurcations of domain walls 

To put this concept in concrete terms, we focus on solutions which are similar to the Bloch 
wall solution, i.e. inversion symmetrical solutions with Di = D2 = Li belonging to the 
fixed points (0,0) and (0,7r) of Li. Applying the above considerations we identify the so- 
lutions by intersection points of or f/(VF]^^) and symmetry lines. 
Symmetric domain wall solutions can be divided into two classes [Qj: Configurations with 
the inversion symmetry point on a lattice site are described as central-spin solutions. They 
correspond to central sequences of B = U oV (p^. 
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The remaining (central-bond) configurations are characterized by an inversion symmetry 
point in the middle between two lattice sites. They correspond to B = U o Li o V (|l^) or 
to B = UoL2oV {^) 

In the following, selected solutions described by central sequences B = Li, L\ , L2 oi L2 
are studied. 

Typical solutions for moderate parameter values = jj = 0.28) can be inferred from fig- 
ure 1 which shows the sector [0, 1] x [0, tt] of the phase space. The sector can be continued 
with the period tt by the symmetry M'^ , and by M~ it is inversion symmetrical. This 
corresponds to clockwise and anticlockwise solutions. 

An intersection point of W^_^ and VF£^ is mapped to the next intersection point but one. 
Therefore two solutions arise from intersection points of W}^ and . Corresponding to 
(21), a heteroclinic set with B = L\ exists which contains the intersection point s of Wf^^, 
Wf^^ and the symmetry line (f) = ^ oi . This solution contains a central spin with </> = f 
(figure 2). 

Corresponding to (p5|) the heteroclinic set with B = Li contains an intersection point b of 
^Li' ^^'^ symmetry hue (p = ^+^''^^'^™ ^ of Hi~^ . The center of this solution is in 
the middle between two lattice sites. (0, 0) is an intersection point of VF^^ and Fix{H2~^); 
so there exists a sudden kink (pi = for i < 0, (pi = ir for i > 0. This solution of type 
B = L2 corresponds to the Bloch wall (^) for j 00. 

A differently shaped domain wall of type B = L2 corresponds to the intersection point PSl 
of L2'^{W}^J and W[^. The predecessor 51 of PSl is an intersection point of L2{Wl_^) and 
the symmetry line (p = ^ of GJ"*", so the centre of this solution is (p = ^ (figure 3). 
If a localized solution arises from one and the same domain on its right and its left side, it 
is called hump. We pick out three simple hump solutions represented by homoclinic orbits 
of the dynamical system. The intersection point C of L2{W]^_^), W^^ and the symmetry 
line ^-^''^^"^■^ of H2~^ indicates a homoclinic B = L2 solution, which connects (0, 0) to itself 
and is centred at (p = tt. 

There exists also a homoclinic solution of B = L2 type which can be seen from the intersec- 
tion point PS2 of L2^{Wf^^) and Wf^^. The predecessor 5*2 of ^5*2 is an intersection point 
of L2(W^") and the symmetry line (p = tt oi G^^; the solution is centred at = vr. The 
origin (0, 0) is also an intersection point of these two lines which corresponds to a trivial 
hump (pi = 0, 00 = ^1 'Pj = 0, i < 0, J > where only one spin is folded down in relation to 
the ferromagnetic state. 

In the following, bifurcations of these solutions will be briefly described. To determine the 
stability of solutions in question the equations of motion for a number (usually 30-50) of 
spins of the domain wall are linearized while the boundaries are fixed. Eigenvalues of this 
matrix are computed numerically. Due to the rotation symmetry of the problem, two eigen- 
values are zero. The solution corresponds to a local energy minimum if all other eigenvalues 
are imaginary. This fixed point becomes attractive if a phenomenological Landau-Lifshitz 
damping term |12| is added. 

The solution corresponding to intersection point b is stable. It exists up to = 0.53, where 
the maximal / of the Li solution (point b) reaches 1. Li and L2 have the same effect at this 
point. It follows that a B = Li solution is identical with & B = L2 solution. Indeed, with 
this parameter value a branch of L2(W}^_^) touches VF£^. So the B = Li solution merges 
into a B = L2 solution with an opening angle larger than 5. This solution is stable as well. 
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For increasing k, the intersection point of L2{W^ ) and is moving down to I = 0,(j) = tt. 
At this point, the clockwise L2 solution and its anticlockwise counterpart merge into the 
sudden kink solution in a pitchfork-bifurcation. In this bifurcation the sudden kink solution 
corresponding to the intersection point (0,0) of Wj^^ and Fix{H2~^) becomes stable. We 
can compute the parameter value of this bifurcation analytically. The slopes of VF£^ and 
of L2{Wj^J at (0, vr) must coincide, i.e. 

-k-V2kTk^= k + V2k+^ 

i-k-V2kT¥ 

must be satisfied, which happens at A: = |. This is the analytical confirmation of the 
numerically determined value in 

The solution containing a central spin s is stable and has a slightly higher energy. It exists 
for all parameter values. For j ^ 00 we get 0i — > for i < 0, (pi = ^ for i = <j)i ^ ir for 
i > 0. 

The solution corresponding to the intersection point S2 appears in a pitchfork-bifurcation 
out of a trivial solution 4>i = 0,(f)Q = tt, = 0, i < 0, j > 0. At the critical parameter value 
the gradients of W[_^ and L2^(W^l^) must coincide. The value is computed analogously to 
(|2^ , we get k = J. For k < 1.28 the trivial hump is unstable, for larger k it is stable. All 
other solutions of figure 3 are unstable. 

Conclusion 

In contrast to the related continuous system, the discrete spin chain has a large variety 
of stationary planar solutions. As a peculiarity of spin systems, a virtual sequence of two 
mappings generates a stationary solution. Consequently a solution is defined both by its 
initial conditions and a generating symbol sequence. Besides a variety of solutions with no 
correlate in the continuous model, the discrete spin chain has several solutions similar to 
the Bloch wall solution. The latter leads to degeneration if the continuum limit is a good 
approximation (weak anisotropy) : Two stable solutions of the discrete system converge to 
the same continuum limit. On the other hand, for strong anisotropy, one of these solutions 
bifurcates into an abrupt kink of 180 degrees. Among the immense number of solutions of a 
nonlinear mapping, chaotic orbits can be identified with stationary states of the spin chain. 
For moderate values of k, some KAM-tori have already vanished but others still ensure 
that the trajectory does not leave the definition area. Because of energetical instability, 
the physical relevance of these solutions should be assessed with care. Even stable solutions 
may be unobservable due to their high energy, but also unstable stationary solutions may 
be relevant in dynamics. Firstly, unstable stationary and periodic solutions contribute to 
the structure of chaotic attractors and are relevant for their construction [l^. Secondly it 



was suggested that unstable stationary solutions influence transient dynamics in the 



process of pattern formation. 
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Figure captions 



Figure 1 

The sector [0, 1] x [0, tt] of the phase space for k = 0.28. It shows W^, L2(M^"), L2^(W^^) 
and all symmetry lines. Here L2(VF") is the image of of the segment (f) G [— tt, 0]. There 
are a central-bond set of intersection points of W'"" and including the intersection point 
6 of H^", H^^ Fix{Hl+) and a central -spin set including the intersection point s of W^, , 
Fix{G'^^). The intersection point B of L2(M^"), and Fix{H2'^) belongs to a central- 
bond solution. There are two other central-spin solutions: One of them corresponds to 
the intersection point PSl, its predecessor SI is a intersection point of L2(W'"') and the 
symmetry line 4> = tt. The predecessor 5*2 of the point ^5*2 is an intersection point of 
L2{W'^) and ^ = 0. The origin of this solution is a pitchfork bifurcation at A; = 0.25 where 
the gradients of and L2^{W'^) at (0,7r) coincide. 

Figure 2 

The angle 4> of three basic domain wall solutions as a function of the spin index N at 
k = 0.28. The kinks corresponding to symmetry lines C^"*" and H'^''^ are energetically 
stable, these solutions are generated only by Li. The sudden kink corresponding to the 
symmetry line G^~^ is caused by a central application oi B = L2. It becomes stable at 
A; = I (weak coupling and virtually isolated spins). 

Figure 3 

The domains ^ = and (p = tt are linked by a solution corresponding to the symmetry 
line G"^"^ (note mod27r). Discernible from opening angles larger than ^, this solution 
contains a central sequence B = L2^. So PSl is generated from its predecessor 51 by L2. 
Similarly, the hump G^^ contains B = L2^, while the hump contains B = L2. At the 
trivial hump G^"*" , only one spin 52 is folded down in relation to the ferromagnetic domain. 
This solution becomes stable for weak coupling {k = 1.28) while all other solutions in this 
picture are unstable. 
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